-
Notifications
You must be signed in to change notification settings - Fork 2
/
sqrt.c
62 lines (45 loc) · 1.35 KB
/
sqrt.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
/*++
toro C Library
https://github.com/KilianKegel/toro-C-Library#toro-c-library-formerly-known-as-torito-c-library
Copyright (c) 2017-2024, Kilian Kegel. All rights reserved.
SPDX-License-Identifier: GNU General Public License v3.0
Module Name:
sqrt.c
Abstract:
Implementation of the Standard C function.
Calculates the square root of a floating-point value.
Author:
Kilian Kegel
--*/
#include <CdeServices.h>
extern double __cdecl __cde80387FSQRT(double x);
/**
Synopsis
#include <math.h>
double modf(double x, double* intptr);
Description
https://learn.microsoft.com/en-us/cpp/c-runtime-library/reference/sqrt-sqrtf-sqrtl
Parameters
https://learn.microsoft.com/en-us/cpp/c-runtime-library/reference/sqrt-sqrtf-sqrtl#parameters
Returns
https://learn.microsoft.com/en-us/cpp/c-runtime-library/reference/sqrt-sqrtf-sqrtl#return-value
**/
double __cdecl sqrt(double d)
{
CDEDOUBLE* pdbl = (void*)&d;
uint64_t di = 0xFFF8000000000000ULL;// 0x7FF8002000000000LL;//DOMAIN ERROR
double* pd = (void*)&di;
do
{
if (0 == pdbl->member.exp) // 0.0?
{
*pd = d; // return 0.0
break;
}
if (0 == pdbl->member.sign) // non-neg?
{
*pd = __cde80387FSQRT(d);
}
} while (0);
return *pd;
}